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ABSTRACT 



Atmospheric turbulence introduces random phase distortions in optical 
imaging systems. The development of new laser and imaging systems requires 
information on the spatial and temporal distribution of this atmospheric 
turbulence. Measurements of the image spread and the jitter induced by the 
atmosphere on an optical system provide two techniques to quantify these 
phenomena. This thesis evaluates a Spectra Sources Lynxx PC Plus charge 
coupled device (CCD) array as an atmospheric turbulence sensor. Data 
acquisition and processing programs were written to measure the image spread 
of a point source and centroid jitter of a point source imaged through the 
atmosphere. Since atmospheric jitter measurements require high image frame 
rates on the order of 200 images per second, a large portion of this thesis 
involved measurements of the times for the CCD detector, interface board and 
IBM compatible computer to perform their tasks. Recommendations for higher 
performance are presented. 
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I INTRODUCTION 



This thesis considers a method used to measure atmospheric turbulence 
using a silicon charge coupled device (CCD) camera. Light from a distant star 
consists of almost perfect plane waves before entering the earth’s atmosphere. 
When the waves traverse the atmosphere, local variations in the index of 
refraction induce warping of the wave fronts. Therefore, a telescope pointed at 
a star will produce a distorted image, whose spatial and temporal characteristics 
depend upon the properties of the atmospheric path followed by the starlight. 
Typically, the atmosphere introduces both image spread and temporal jitter, so 
that a CCD array placed at the telescope’s focal plane will reveal a highly 
unstable and complex irradiance pattern. It has been found that both a high 
frame rate (> 200 per second) and correspondingly short exposure time (< 5 
ms) are necessary to avoid under-sampling the dynamic variations of typical 
images. 

The goal of this thesis was to find a way to measure atmospheric turbulence 
with a CCD camera by optimizing the speed and efficiency of the camera 
system’s software and hardware. The software was written in the C language. 
Software components that can affect system speed include both the various 
algorithms used and the overall efficiency of the compiled program code. Three 
compilers’ outputs were compared for executable code speed and utility: 
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Microsoft Quick C 2.5, Borland C++ 2.0 and Borland Turbo C++ 1.0. Fixed 
hardware components included a Texas Instruments TC211 CCD image sensor 
in a Spectra Source PC Lynxx Plus camera attached to an 8 bit IBM PC 
compatible interface card. The camera system was installed in two IBM 
compatible personal computers, one using an Intel 80386 and the other a 80486 
processor, and the overall speeds were compared. Finally, the camera-computer 
interface and the CCD image sensor were analyzed to determine whether such 
a low-cost commercial system would be suitable for turbulence measurements, 
or whether a customized device would have to be fabricated. 
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II BACKGROUND 



A. ATMOSPHERIC TURBULENCE 

Turbulence in the atmosphere causes many problems. In astronomy, it 
introduces distortion that obscures image detail. Since exoatmospheric use is 
not practical for most telescopes, knowing the spatial and temporal distribution 
of turbulence is essential in the planning stages of a new facility or in diagnostic 
evaluation of test systems. Figure 1 sketches a wave front incident on a turbulent 
region [Ref. 1]. Random phase fluctuations in the index of refraction field 
produce scintillations and image blurring. Twinkling or scintillation arises from 
the interference of starlight that traverses multiple paths through the atmosphere. 
Image blurring arises from the reduction in spatial coherence of the wave front. 



B. PARAMETERS FOR THE MEASUREMENT OF TURBULENCE 

The atmospheric index of refraction depends on both temperature and 
pressure. Since pressure fluctuations disperse rapidly, at the speed of sound, 
temperature fluctuations are the main contributor to atmospheric optical 
turbulence. Structure functions provide a way to quantify the statistics of 
atmospheric turbulence. For propagation through the atmosphere, optical 
parameters will involve an integral of the refractive structure function along the 
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Figure 1 . Optical atmospheric turbulence distorts a coherent phase front, 
optical path. The primary optical parameter is the spatial coherence length, r„. 
It represents the transverse autocorrelation length of the electromagnetic field. 
A related parameter is the isoplanatic angle, 0^. It represents the angle between 
two different paths where the Strehl ratio of an ideal adaptive optics system is 
within e ’ of perfect correction [Ref. 2]. Another parameter, the Greenwood 
frequency, f^, represents the electrical bandwidth needed to remove atmospheric 
phase distortions with an adaptive optical system. 
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1 . structure Function 



Atmospheric turbulence produces localized variations in the index of 
refraction along an optical path. Tatarski [Ref. 3] defines the structure function 
by 

= ([xfFj) -x(Fj^) ] ^>, (1) 

where < > represents an ensemble spatial average, T, and ~r^ represent the 
location of two points in space, and x represents an atmospheric parameter of 
interest, such as temperature or index of refraction. Assuming the turbulence is 
isotropic, homogeneous and incompressible [Ref. 4] the Kolmogorov theory of 
turbulence [Ref. 5] shows that 

D^=clrli\ ( 2 ) 

where C/ is a structure parameter and r,g is the distance between the two points 
in space. By combining equations (1) and (2), the thermal structure parameter, 
which characterizes the mean squared temperature difference between two 
points in space is [Ref. 3], 



rii 

where Tg-T, is the temperature difference and r,j is the distance between the two 

« 

points in space. There is a similar expression for the index of refraction structure 
parameter, Taking the partial derivatives of the atmospheric index of 

refraction [Ref. 3], 
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n-1 



79x10'^ P 
T 



(4) 



where P is the atmospheric pressure in mbar and T is the Kelvin temperature, 
and assuming isobaric turbulence allows us to write in terms of C/ [Ref. 6], 

= ( 79 X10-^P ^2 

The pressure is assumed to be constant, since small random pressure 
differences disperse rapidly. 

2. Spatial Coherence Length 

The spatial coherence length, r^,, as stated above, measures the 
transverse autocorrelation length of a wave front. It represents the effective 
aperture diameter of a diffraction-limited optical system with a similar angular 
resolution as the system under study. Typical values of r^, range from two or 
three centimeters for high turbulence to as much as thirty centimeters for low 
turbulence, when measured from the earth’s surface upward. When observing 
a star from the ground with an optical device, the component of r^ is 
cumulative along the optical path, even though most of the degradation of the 
wave front occurs close to the ground. Integrating the optical turbulence Cn^(z) 
along the path [Ref. 6], 

^ k^sec4>j^^C^(z) W{z) dz] , (6) 

where k is the wave number(27r/A.), (p is the zenith angle, L is the vertical path 
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length and W(z) is a weighting function. Typically W(z) has one of three forms: 
for a plane wave, W=1; for a spherical diverging wave, W=(z/L)®'®; and for a 
converging wave, W=(1-z/L)®'^. The factor of sec 0 compensates for the slant 
path through a horizontally stratified atmosphere. 

Temperature fluctuations along a vertical path have been measured by 
two means: by microthermal probes carried by a meteorological balloon and by 
an echo sounder [Ref. 7 and 8]. Both of these methods provide credible values 
of r^,. However, there is a more direct method for calculating the coherence 
length at the earth’s surface using an optical point source irradiance distribution 
[Ref. 9]. The discrete irradiance values from a two dimensional irradiance 
distribution represent the point spread function, P(x,y). The line spread function 
for X, L(x), is the summation of all of the y irradiance values from P(x,y) for each 
X. L(y) is a similar summation of the x values for each y. Then take the Fourier 
transform of L(x) or L(y) to get the optical transfer function (OTF). Assuming the 
turbulence is isotropic, that is symmetric with respect to rotation about the path, 
and that it is laterally stationary, the absolute value of the OTF equals the 
modulation transfer function (MTF). For an optical system imaging a star the 
observed quantity, MTF^, is a product of the MTF’s of the source or star, the 
atmosphere and the instrument [Ref. .9], 

MTF^ = MTF^'MTF^-MTF^. ( 7 ) 

MTFj for a star corresponds to a point source at infinity and is unity. MTFj is a 
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measured instrumental transfer function. To determine the atmospheric 
component, MTF^, from the observed quantity, MTF„, we divide MTF^, by the 
measured instrumental transfer function, MTFj. 

Fried shows that atmospheric MTF^ is a function of the atmospheric 
wave structure function, which is an integral of the index of refraction structure 
parameter along the path [Ref 2]. MIF^ reduces to a form [Ref. 9], 

- 3.44 

MTF^ = e , (8) 

where R is the effective focal length, A. is the wavelength of light and v is the 
spatial frequency. To determine the atmospheric parameter r^, we took a one 
dimensional fast Fourier transform of a line spread function of a star image. 
Dividing this by the known instrument function gives the spatial spectrum similar 
to Equation (8). The next step was to determine the e ’ frequency for this 
experimental spatial spectrum from which [Ref. 9], 

= 2.1 A. i? v^, (9) 

where was the e ’ frequency of the MTF^. 

Yet another way to calculate r^ uses image centroid jitter, also called 
beam wander, by calculating the x and y centroid standard deviation from a 
sequence of point source images. It is equal to aF, where a is the atmospheric 
angle of arrival fluctuations and F is the optical system’s effective focal length. 
The angle of arrival is sketched in Figure 2. Mathematically, angle of arrival 
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Figure 2. The angle of arrival is the least-squares plane that fits the electric field 
vector. 

fluctuations are equivalent to least-squares planes that fit the electric field across 
an aperture. In terms of the index of refraction structure parameter [Ref. 5 and 



where D is the diameter of the aperture and L is the optical path. Combining 
equation (10) with equation (6) for a plane wave, 



To calculate r„ from the centroid motion requires a very high image frame rate to 
"freeze" the atmosphere. Characteristically at least 200 frames per second are 
required in order to avoid aliasing that would otherwise underestimate the true 
centroid variance. 



10 ], 




(10) 



= [1 . 418 sec([) . (11) 
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3. Isoplanatic Angle 

The isoplanatic angle measures the angular coherence in the vicinity 
of an object. It defines a cone that constrains the angles over which an adaptive 
optical system will provide valid correction. In terms of [Ref. 1], 

00 = [2 . 91 f (z) dz] , (12) 

where the parameters are the same as those of r^,. Simple instruments exist to 
measure 0^, [Ref. 11]. For this reason, only the coherence length is addressed 
in this thesis. 

4. Greenwood Frequency 

The Greenwood frequency was introduced in an earlier paper by Darryl 
Greenwood [Ref. 12 and 13]. It is a mean temporal frequency of an adaptive 
optics system and depends on the wind velocity, V, and along the path. For 
a Kolmogorov turbulence spectrum, the Greenwood frequency is [Ref. 13], 

f = dz]^/^ (13) 

where X is the wavelength of light, z is the distance along the propagation path, 
V(z) is the wind speed and L is the distance from the source to the receiving 
optical device. 
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C. CHARGE COUPLED DEVICE (CCD) 



An atmospheric coherence length sensor needs an imaging detector to 

» ..j 

measure the point spread function. The speed, reliability and availability of 
charge coupled devices (CCD’s) makes them worthy of consideration. The basis 
of a CCD is a metal-oxide-semiconductor (MOS) capacitor forming potential wells 
and channels that comprise light-sensitive pixels and read-out registers. Figure 
3 is a sketch of the MOS capacitor structure [Ref. 14], Electrons excited by the 
photoconductive effect are trapped in potential wells formed under the positively 




Figure 3. The metal-oxide-semiconductor (MOS) capacitor structure, the basis of 
the charge coupled device (CCD). 

charged metal gate contact. The number of electrons trapped is proportional to 
the integrated irradiance during an exposure. 
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There are several methods to read-out the device. Voltages with different 
phases placed on the gates transfer charge from pixel to pixel or potential well 
to potential well. A set of electrodes connected together is called a phase. 

I 

Figure 4 shows an example of a three phase device with three pairs of gates (G) 
and three lines (L) connecting the electrodes [Ref. 14]. A positive electrical pulse 

L, 

M 

0 

(A 

e 
8 

I: 

1 

Figure 4. Three pairs of gates (G) are connected with the 
lines (L) to form three different phases. 

from a clock forces a transfer of charge from one phase of pixels to the next. 
Clocking a phase of cells reduces the barrier between pixel columns allowing 
charge transfer. A bucket of water running down several steps with boards or 
barriers on each step is a visual representation of this process. The water flows 
down to a lower potential each time the board is lifted up. The number of 
phases used to transfer out charge packets represents the number of steps. 

There are two methods used to insure that efficient charge transfer occurs 
in one direction. The first uses multiple phases to separate the charge packets. 
The Texas Instruments TC21 1 [Ref. 15] uses the second method which is to force 
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an asymmetry in each well using an ion implantation zone between each pixel. 
The clocking method is the monophase mode, also referred to as a 1 + V 2 phase 

•o 

mode, sketched in Figure 5 [Ref. 16]. One phase has an intermediate voltage 

I 

level while the voltages applied to the other phase vary on both sides of this 
level. 

The charge packets transfer down channels. There are two types of 
channels, surface and buried. The water analogy also works here. In a buried 
channel, there is less loss of charges, similar to a pipe. A surface channel, 
however, has a higher capacity, like a canal, but this CCD has higher noise from 
surface imperfections. 



01 



o? 



01 







Figure 5. The monophase clocking method holds a phase at an intermediate 
level while the voltages are applied on either side. 
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During the read-out of a CCD, image smear will occur since light will still 
produce charge in the moving CCD rows, unless a shutter is incorporated. 
Different techniques exist to reduce image smear during read-out. Frame transfer 
devices overcome this disadvantage by effectively having two arrays, one for 
image exposure and a second for storage. Moving charges from the image area 
to an opaque, masked set of pixels quickly allows another exposure to 
commence while reading the storage area. Figure 6 shows two examples [Ref. 
14 ]. 
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Figure 6. Two schemes for CCD read-out. Top, the interline method; bottom, 
frame transfer method. 
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Ill DISCUSSION 



A. EQUIPMENT 

This thesis investigated a Spectra Source Lynxx PC Plus CCD Imaging 
System [Ref. 17 and 18]. The original intent was to use this CCD camera with 
a telescope and personal computer to develop software to measure the 
atmospheric coherence length, r,,. During this process it became clear that the 
off-the-shelf CCD system could not provide sufficient frame rates to achieve the 
intended goal. Consequently, the task shifted to determine the causes of the 
sluggish acquisition rates. The CCD detector array, computer interface card and 
the IBM compatible PC were investigated. 

1. CCD Camera 

The camera contains a 192 X 165 pixel Texas Instruments TC21 1 CCD 
image sensor. Its specifications are in Appendix A [Ref. 19]. Figure 7 is a 
photograph of the sensor centered inside the camera. The other components 
inside the camera head translate the clock signals from the interface board and 
amplify the CCD output. 

Figure 8 shows a functional block diagram of a TC21 1 CCD image 
sensor [Ref. 19]. Key components include the silicon matrix of pixels and a set 
of gates to shift the charge. The image area gate (lAG) performs a parallel shift 
of all rows for each clock pulse. As a row enters the output serial register the 
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Figure 7. Photograph of inside the CCD camera head of the Spectra Sources 
Lynxx PC+ system. 

serial register gate (SRG) shifts the row to the charge detection amplifier. The 
TC211 includes an antiblooming gate (ABG) that should not be used since it 
introduces a severe nonlinear photo response. 

The method used to acquire and to read-out an exposure determines 
the CCD sensor’s read-out rate. The electrons collect in potential wells, or pixels 
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Figure 8. Lynxx PC+ image sensor CCD, Texas Instruments TC211 block 
diagram. 

after exposing the image area to light. The charges then shift down to the output 
register as rows, through 165 timing cycles. The output register then shifts the 
charges 210 times for each row. This cycle of 210 shifts includes 12 dark pixels 
for a dark reference and 6 dummy pixels used to transfer the charges out of the 
register. Figure 9 shows the charge transfer process [Ref. 19]. Another factor 
to consider for an image sensor chip is the noise. The noise equivalent signal 
for the TC211 is 150 electrons. It depends on (kTC)’'® /q, where C is the 
capacitance of the read-out charge collector, k is the Boltzmann constant, T is 
the absolute temperature and q is the charge of an electron. 

The integration time needed to achieve a particular signal-to-noise ratio 
depends on both the telescope optics and the object’s radiance. It will be a 
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Figure 9. TC21 1 image area transfer process, uses a virtual phase. 

constant time added to the read-out time regardless of the number of pixels used 
or transfer times. The image area gate (lAG) and serial register gate (SRG) 

rate. 



maximum pulse rates determine the sensor read-out 



The maximum rate for the lAG parallel row transfer gate is 1 .5 MHz. 
This corresponds to 0.667 microseconds (^s) per row. Clocking the lAG shifts 
the whole array down by one row. This can be used to clear the CCD. To clear 
all 165 rows requires only 110 \is. Since an lAG clock pulse clears the serial 
register, the lAG should not be clocked any faster than the SRG can read-out a 
row of pixels. 

The maximum rate for the SRG single row shift gate is 10 MHz. This 
corresponds to 0.1 ^ls per pixel. Figure 8 shows the pixel read-out order. Six 
dummy pixels will appear before any valid data and the twelve trailing dark pixels 
can be ignored. At the maximum clock rates, the ideal minimum time required 
to read an entire row of pixels is 21 0 pixels X 0.1 ^lS = 21 |iS. The minimum time 
needed for all 165 rows is 165 rows X 21 |is = 3465 ps. 

The total read-out time must include the serial shift time and two 

parallel array shifts. One set of lAG parallel shifts are needed to move each row 

into the output shift register and a second set must clear the CCD before the next 

exposure. Since the output register cannot be read during a row transfer (lAG 

clocking), using all the pixels, the ideal minimum read-out time is 3.685 ms. 

Placing the image in the lower left portion of the array can increase the frame rate 

since only a portion of the array needs to be read-out. After shifting the rows 

« 

with information to the output register, the lAG clocking can increase to its 
maximum rate to clear the CCD and to commence another integration. The six 
leading dummy pixels will always add some time to each serial read-out time, but 
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the twelve trailing dark pixels could be chopped off. As an example, a 60 X 60 
image would take 2 X 60 rows X 0.667 \is + 66 pixels X .1 ^s X 60 rows for the 
serial read-out for a total of 476 ^s. There are 66 pixels per row because of the 
six dummy pixels. 

As mentioned earlier, the integration or exposure time adds to the read- 
out time to produce the total measurement time. The total time could be reduced 
to only the integration time if the CCD were read-out while integrating the next 
set of data. This reduces the minimum total time to the longer of these two 
processes. Frame transfer devices as shown in Figure 6 have a complete extra 
set of pixels to hold the data for read-out. The interline method is not practical 
for a star image because the gaps at the transfer pixel rows reduce the exposed 
area by a factor of two. Each pixel would lose half of its photons. This is 
unacceptable for this application. Masking off half of the array is a possible 
solution, it would not work with the TC21 1 chip. Clocking the lAG parallel row 
shift register would produce streaking in the upper half exposed pixels while 
reading the masked pixels. 

2. Interface Card 

The Spectra Source computer interface card contains the circuitry 
needed to transfer information from the CCD array to the computer bus. It does 
an eight bit parallel transfer to an IBM compatible computer. Figure 10 is a block 
diagram of the components of this card. The CCD sensor was discussed in 
Section 1 . The sample-and-hold (S&H) amplifier chip converts video voltage from 
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1 Pixel 




Figure 10. Spectra Sources Lynxx PC+ interface card block diagram. 

Section 1. The sample-and-hold (S&H) amplifier chip converts video voltage 
from one pixel into a voltage pulse. It is a NEC SE/NE5537 [Ref. 20]. The 
minimum acquisition time is less than 4 ^lS. It is controlled by pulses from the 
chip below it, a 74LS123P [Ref. 21], which is a dual retriggerable one-shot. Its 
pulse switching characteristics are in the nanosecond (ns) range. A 1200 pF 
capacitor and 12 k resistor determine the present pulse length of 5180 ns. The 
control chip characteristics did not limit the speed of the S&H chip. 

The analog digital converter (ADC) chip converts the voltage to a 
positive 1 2 bit integer, which has a maximum value of 4095. This system has a 
Maxim MX7572 ADC chip [Ref. 22]. Its minimum conversion time is 5 ^is. It, as 
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and-hold and digital conversion time is 9 jis for each pixel. For a 60 X 60 pixel 
image, the total S&H and ADC time is 32.4 ms. 

3. Computers 

Two different computers were used to investigate the dependence of 
the frame rate on computer characteristics. The first was a Compaq 386, 20 
megahertz machine. It had a 16 megabyte RAM and a 60 megabyte hard disk. 
The second was a Dell 486, 33 megahertz machine. It had a 16 megabyte RAM 
with a 350 megabyte hard disk. 

B. SOFTWARE 

The Lynxx system came with software programs designed for amateur 
astronomical observations [Ref. 17]. This thesis used a separate set of data 
acquisition program modules written for NFS by Spectra Sources [Ref. 18]. It 
was necessary to modify these modules during this thesis research. The 
software was written in the C language. Microsoft Quick C 2.5, Borland Turbo 
C + + 1.0 and Borland C++ 2.0 compilers produced executable code that were 
tested to compare their relative speeds while collecting and processing data. 
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IV RESULTS 



The first task involved writing the modulation transfer function, centroid jitter 
and display programs. After it become clear that the Lynxx PC + CCD system 
was too slow for effective centroid measurements, the second task developed 
programs to analyze the Lynxx image acquisitions times. The data acquisition 
functions were timed using different computers and code from different 
compilers. The manufacturer’s maximum speed specifications for both the CCD 
camera and the interface card were compared to the actual measured times. 
After analyzing these results, replacements were recommended. 

A. PROGRAMS 

A modulation transfer function program and a jitter program were written to 
compute the atmospheric coherence length, r^. An artificial star program was 
developed to test the MTF program. The following sections describe these 
programs. 

1. Modulation Transfer Function 

This program calculated the coherence length, r^,, from a stellar image 
intensity distribution. It subtracted the background from the star image, 
calculated the line spread function, found the centroid and calculated the fast 
Fourier transform. It then calculated r,, using equation (9). The calculations were 
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done in the x and y dimensions separately, giving a quality check, since actual 
image turbulence is nearly always isotropic. With the image centroid data the 
telescope could operate in an auto track mode, a feature that comes with the 
Lynxx software package. The program plotted the modulation transfer function 
as a visual check. The program separates into modules for ease of future 
revisions. This program is in Appendix B. 

2. Artificial Star Program 

This program created a perfect exponential star image for the MTF 
program. This provided a convenient test of the MTF program, since the FFT of 
an exponential is also an exponential. This program is in Appendix C. 

3. Jitter 

This program calculated r^, from the centroid jitter. It first removed the 
background by subtraction. After this, it computed the image centroid jitter using 
the standard deviation of the image, according to equation (11). This program 
also calculated the line spread functions and plotted them on the screen so that 
proper telescope tracking was verified. Appendix D contains this jitter program. 

B. TIME MEASUREMENTS 

Programs were written or modified to compare frame rates attainable for 
different subframe sizes. The computer clock measured the time differences. 
Since the MS-DOS system clock had a 55 ms resolution, many of the times were 
done for 100 cycles of the function of interest. All the times used were at least 
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one second long, resulting in an accuracy within 1%. An array of 5 or more trials 
ensured consistency. Appendix E is an example of the code used for measuring 
the times. 

1. Computer 

We found that the type and speed of the computer did not contribute 
significantly to the camera frame rate. Measurement time included the clearing 
of the CCD, the integration, and digitization times. With an image size of 60 X 
60, there was less than V 2 of a frame per second difference between 386 and 
486-based computers. Figure 1 1 shows a comparison of frame rates verses 
subframe size for the two machines. Including calculations for the line spread 
functions and backgrounds, the frame rate improvement was about 3V2 frames 
per second for a 60 X 60 image with a 486 computer over that achievable with 
a 386 machine. This is shown in Figure 12. The mathematical calculations were 
timed to determine their contributions to the total frame rate. These calculations 
could have been done separately if they had slowed the frame rate, but they 
were actually performed in real time. 

2. Compiler 

It was found that the choice of compiler did not significantly influence 
the speed of the measurements. Figure 13 shows there was no difference in the 
frame rate for the code produced by three different compilers to acquire data. 
Figure 14 shows that after adding the mathematical calculations to the acquisition 
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COMPARISON OF PROCESSORS WITHOUT CALCULATIONS 
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Figure 11, Spectra Sources Lynxx PC+ frame rates using IBM PC 
386 and 486 computers. Times include clear, integration and 
digitization, using code compiled under Turbo C++. 
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COMPARISON OF PROCESSORS WITH CALCULATIONS 
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Figure 12. Spectra Sources Lynxx PC+ frame rates using IBM PC 
386 and 486 computers. Times include clear, integration, 
digitization and mathematical calculations, using code 
compiled under Quick C. 
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COMPARISON OF COMPILERS WITHOUT CALCULATIONS 
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Figure 13. Spectra Sources Lynxx PC+ frame rates using Turbo 
C++, Borland C++ and Quick C compilers. Includes clear, 
integration and digitization times, using a Dell 486, 33 MHz 
PC. 
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COMPARISON OF COMPILERS WITH CALCULATIONS 
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Figure 14. Spectra Sources Lynxx PC+ frame rates using Turbo 
C++ / Borland C++ and Quick C compilers. Includes clear, 
integration, digitization and mathematical calculation times, 
using a Dell 486, 33 MHz PC. 
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time, the code produced by Quick C was about % of a frame/second slower than 
the code from Borland C + + or Turbo C++ for a 60 X 60 image size. 

3. Functions 



The times for each measurement were divided into four steps. Table 
1 shows these steps. 

TABLE 1 



Step 


Image size 


Time 


Cl©3r CCD 


Must clear entire array 


constant 6.45 ms 


Integration 


Independent of size 


used 5.0 ms 


Digitization 


60 X 60 


86.1 ms 


Calculations 


60 X 60 


40.2 ms 


Total 


60 X 60 


137.75 ms 



The digitization and calculation times depended on the size of the image, while 
the clear and integration times were constant for each measurement. This is 
shown in Figure 15. Clear and digitization times depended on the image sensor 
and interface card. 

C. HARDWARE RATE COMPARISONS 

This section compared the maximum design speeds of the 
chips in both the camera head and interface card to those actually measured. 
A general time reference was 200 frames/second which corresponded to 5 ms 
per frame. 



31 



TIME VS FRAME SIZE 
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Figure 15. Spectra Sources Lynxx PC+ function times for a 
Dell 486, 33 MHz PC. Clear and integration times were 
constant at 6.45 and 5.0 ms respectively. Digitization and 
calculation times depended on the subframe size. 
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1. CCD Camera 



The measured time for the clear function was 6.45 ms. The present 
sensor, TC211, is capable of clearing the entire array in 110 ^lS. All rows need 
to be cleared even when reading a sub array such as a 60 X 60 image. This 
function was much slower than the sensor capability. It is not known why the 
clear times are so slow. 

The present sensor could read a 60 X 60 image in 436 jis at maximum 
design speeds. The measured digitization time, which included the read-out 
time, sample-and-hold (S&H) and the digital conversion (ADC) times was 86.1 
ms. Since this was much larger than 436 ^s, the acquisition time depended 
primarily on the S&H and ADC processes. 

2. Interface Card 

The interface card controls the S&H and ADC times. Their maximum 
speed specifications would result in processing a 60 X 60 image in 32.4 ms. This 
is a significant portion of the observed 86.1 ms digitization time. We therefore 
concluded that these two chips were the primary contributors to the frame rate 
and were not operating at their maximum speed specifications. The sensor itself 
can hypothetically operate at significantly higher speeds. 
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D. RECOMMENDATIONS 



1. CCD Camera 

To increase the frame rate requires an image $ensor that uses a frame 
transfer method with separate image and storage arrays. The Texas Instruments 
TC277, a 735 X 580 pixel CCD image sensor is a proposed replacement sensor 
[Ref. 19], Figure 16 shows its set-up. The image area is 699 X 288 which is more 
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Figure 16. Recommended image sensor CCD, Texas Instruments TC277 block 
diagram. 

than big enough for our applications. The image sensing area and image 
storage area have different gates. The image area gate (lAG) shifts the rows from 
the sensing area to the storage area one row per clock pulse. The storage area 
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gate (SAG) shifts the rows down to the serial register at one row per dock pulse. 
The three serial register gates (SRG) shift out the pixels in each row. 

The storage area gate (SAG) maximum clock rate is 3.34 MHz, which 
corresponds to 0.299 ^xs per row. The serial registers gates (SRG) maximum 
clock rates are 4.46 MHz, which corresponds to 0.224 ^lS per pixel. For a 60 X 
60 image, noting that there are eleven dummy pixels per row, the read-out time 
is 60 rows X 0.299 ^lS + 71 pixels X 0.224 ^iS X 60 rows for a total of 972 ^iS. 
This is done in parallel with the integration. Since this is much less than any 
expected integration time, the integration time and transfer time of data from the 
sensing area to the storage area become the limiting factors for frame rates. 

The image area gate (lAG) maximum clock rate is 3.34 MHz, which 
corresponds to 0.299 jis per row. The dear time would then be 288 rows X 
0.299 ^is or 86 ^lS. This dear time plus the number of rows X 0.299 ^lS plus the 
integration time will be the measurement time for each set of data. For an image 
with 60 rows, the measurement time would be only 86 + 17.9 = 104 ^s, plus the 
integration time. For the TC21 1, the minimum clear and 60 X 60 image read-out 
time is 110 -1- 436 = 546 ^iS. The TC277 is over 5 times faster. 

It should also be noted that the number of pixels in each row of the 
TC277 can be the maximum available (735) without affecting the speed, due to 
the parallel read-out and integration times. The noise equivalent signal is only 
25 electrons, a 6-fold improvement from the TC211’s 150 electrons. Replacing 
the image sensor would require modifications to the interface board. 
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Using such an improved sensor, the bulk of the data could be 
transferred while collecting the next data set. What remains is to get the 
sampling and digital conversion times fast enough to be completed during the 
integration time. 

2. Interface Card 

The first step in the digital conversion is the sample-and-hold (S&H). 
A proposed replacement chip for this is Analog Devices AD781 [Ref. 23]. It 
provides a 700 ns conversion time. This is over 5 times faster than the present 
SE/NE5537 S&H chip’s 4 ^s. For a 60 X 60 image, the AD781 S&H time reduces 
to 2.52 ms from 14.4 ms for the present chip. The present control chip, 
DM74LS123P, for the S&H chip does not need replacement. Changing its 
capacitor to 55 pF and resistor to 25 k would adjust its pulse width to 700 ns 
[Ref. 21]. This is the only revision necessary to replace the S&H chip. 

A proposed replacement chip for the ADC is a Maxim Maxi 62 [Ref. 22]. 
This chip is almost identical to the present one, except that its conversion time 
is 3 \iS instead of 5 jis. For a 60 X 60 image, the improved ADC time would be 
10.8 ms. The Maxi 62 could replace the Mx7552 with little revision to the rest of 
the circuit. 

There are some faster S&H ^nd ADC chips available if the board were 
rebuilt. Datel has a S&H chip with a time of only 25 ns and a ADC chip with a 
conversion time of 350 ns [Ref. 24]. A better option is to have the S&H and ADC 
processes combined in one chip. Of those available, the Datel ADS1 18 [Ref. 24] 
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has a throughput rate of 5.0 Mhz, which corresponds to only 200 ns per pixel. 
The DM74LS123P control chip will still work by using a 60 pF capacitor with a 5 
k resistor [Ref. 21]. Besides providing fast tinning, combining these two chips 
would also simplify the surrounding circuit. For a 60 X 60 image, the S&H and 
ADC time would be 720 ps for the ADS118. 

E. SUMMARY 

Table 2 summarizes all of the relevant times. The components are replaced 
going from left to right across the table. They are in order of the complexity of 
replacement. Replacing the CCD sensor would be the significant step requiring 
major modifications to the interface board. There are also parallel times after the 
CCD sensor is replaced. The processes of clearing, integration and image area 
gate (lAG) transfer occur parallel to the read-out and digital conversion. The 
slowest of these times will limit the frame rate. The S&H and ADC times are the 
limiting factors. For this reason and the necessity to modify the board after 
replacing the CCD sensor, it is best to build a new system. 
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'The lowest of these two rates will limit the system. These two rates reflect the parallel times introduced by the 
frame transfer CCD. 



V CONCLUSIONS AND RECOMMENDATIONS 



A CCD camera can measure atmospheric turbulence by measuring the 
image spread of a point source and by measuring the centroid jitter induced by 
turbulence. The Spectra Sources Lynxx PC+ system evaluated in this thesis can 
measure the coherence length using the point source image spread technique, 
but it is much too slow for jitter measurements. Centroid motion and Greenwood 
frequency measurements need a sample rate of 200 Hz or more to avoid 
undersampling the atmospheric dynamics. Achievement of these rates requires 
a frame transfer CCD, so that parallel image exposure and frame read-out are 
possible. Data processing can be performed during an exposure if the time 
needed is sufficiently short, or after digitizing a series of exposures. The Lynxx 
PC-i- system, which provides 12 frames/second, is much too slow for this. Since 
both the CCD sensor and digital conversion components need replacement, it 
is best to design and build a new system. The recommended CCD image sensor 
is a Texas Instruments TC277 [Ref. 19]. The recommended replacement for the 
sample-and-hold, and analog digital converter chips is a combined sampling 
analog-to-digital converter chip, Datel ADS118 [Ref. 24]. The atmospheric 
coherence length programs developed in this thesis should work with a new 
system. 
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APPENDIX A TEXAS INSTRUMENTS TC211 CCD IMAGE SENSOR 



TC211 

102 X 165 PIXEL CCD IMAGE SENSOR 

D32/V JANUARY 1890 



Full Frama Operalion 

Anlibioomtng Capabilily 

Single Phase Clocking for fforizontal and 
Veilicat Transfers 

Fast Clear Capabilily 

Dynamic Range . . . 60 dB Typical 

High Blue Response 

High Photoresponse Uniformity 

Solid Stale Reiiabfiily Wilh No Image 
Burn in, Residual Imaging, image 
Dlsloilion, image f ag, or Microphonics 

6 Pin Dual In line Ceramic Package 

Square Image Area: 

- 2640 |tm by 2640 pm 

- 192 Pixels (if) by 165 Pixels (V) 

- Each Pixel 13 75 pm (if) by 16 pm (V) 

descrlpllon 

11)6 1C2I1 is a full frame chaige coupled dovice (CCD) image sensor dosiyned specifically ior liidusiriai 
applicalions icrtuiiing luggocfncss and small size 1 Ire image sensing area is conligtired Into 1 65 horizontal lines 
each coniaititng 192 pixels Iwcivo addilional pixels aie provided al the end ol each line lo establish a dark 
releience and line clamp The anlrlilooining loalnre is oclivulcd by su{)plying dock pulses to Ihe aniiblooming 
gale, an integral pail ot each image sensing element The charge is converted to signal voltage at 4 pV per 
electron by a high pertormance siruclnre with biiiit in aiilomatic rosot and a vollage reference generator The 
signal Is liiriher bulleied by a low noise I wo slage soiirco luilowor amplilier lo provide high outpul drive 
capabilily 

1 he 1 C2 1 1 Is supplied in a 6 pin dual In line ceramic package approximaloly 7,5 mm (0 3 In ) square ar>d Is 
chai acterized for operation from - 1 0 X to 45 '*C The glass window can be cleaned using any standard method 
lor cleaning optical assemblies or by wiping the siiilaco wilh a cotton swab soaked In alcohol. 
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TC211 

192 X 165 PIXEL CCD IMAGE SENSOR 



D327I. JANUARY 1990 



functional block diagram 
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pin functional description 
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functional description 

The Image sensing area consists ot 165 horizontal Image lines encti containing 192 ptrotosensltlve elements 
(pixels) Eachpixetls 13 75pm(trorlzonlal)by 16 OOnrn (veitical) As tight enters Hie silicon In the Image sensing 
area, liee elections are generated and collected In potential wells (see Figure 1) During (tils time, the 
^nllbloomlng gate Is activated by applying a burst ol pntses every horizontal blanking Interval. This prevents 
blooming caused by the spilling ol charge from overexposed elements Into nelglrborlog elements. The 
anilbloomtng gate Is typically held at a mid level voltage di irlng readout The qrrantity ol charge collected In each 
pixel Is a linear lirnctlon ol the Incident tight and Hie exposure time. Alter exposure and imder dark conditions, 
the charge packets are transferred from the Image area to the serial register at tlic rate ol one Image Hue per 
each clock pulse applied to tire Image area gate Once an Image line has been transferred Into the serial register, 
the serial register gate can be clocked irntit all ol the charge packets are moved out ol the serial register to tlie 
charge delation node at the ainplilier Inpul 

There are 12 dark pixels to the right otthe 192 linage pixels on eactt Image line Tliesedark pixels are shielded 
from Incident tiglit and the signal derlvf*d from them can be used to generate a daik relerence lor restoration 
ot the video black level on ttie next Image tine 
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TC211 

192 X 165 PIXEL CCD IMAGE SENSOR 



03221, JANUARY 1990 

Each clock pulse a|)plie(i to (tie itnaye area gale causes an anluniattc last clear ol tfie 1G2 image pixels and 12 
dark pixels oi the serial register before the next Image hno Is Iransierred into tiie serial register (Note that tiie 
six dummy pixels at ttie lioitl of tiie serial register, wtiich are used to transport charge packets from tlie serial 
register to tire amplilior Input, are not cleared by tiie Imago area gate clock ) The ^automatic fast clear feature 
can be used to Initialize the Image area by transferring all 165 image lines to tl«e serial register gate tmder dark 
conditions witliout clocking tfie serial rogisler gale 
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Figure 1. Charge Accumulation and Transfer Process 
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Figure 2. Timing Diagram, Noninterlace Mode 



absolute maximum ratings over operating free-aIr temperature range (unless otherwise noted) 



Supply voltage range (ADD) (see Nolo I) 0 V to 15 V 

Clock vollage range (or lAG. SnO. ABG lei mlnals -15Vlo5V 

Operating liee air lempeialuie range -30Xto85'’C 

Stoiage lempeialure range -30 C lo OS'C 

I end temperature 1.0 mm (t/10 Inch) Irom case (or 10 seconds 20OX 

rjOlE I AH vot1n^|« voVtos nro w»Hi ic*»po(1 lo 11*0 nio lofMihml 
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Instruments 
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recommended operating conditions 





MIN NOM MAX 


UNIT 


Sif)pty voHagft. AOB 


It 12 13 


V 


SiAjslialo bias voHaga 


0 


V 


Clock vo11jc)tjl 


Imui\}u etoo gale 


High level 


0 1 2 


V 


liilofinotkiik) lovol^ 


-to ^5 2 


1 ow level 


^10 -9 6 -es 


Soiiai rugiblui gala 


liMjh level 


t 5 2 26 


low level 


-10 -es -e6 


Afild^loofiiing gale 


lligli level 


25 35 4 


liilufiiuKkala level^ 


-3 -2 5 -2 


1 owtovol 


-0 -7 -0 


Clock i<il()S 


lAG 


t 5 


MHz 


SRG 


to 


ABG 


2 


M 


Time inlorval, SRG 1 lo IAGT 


70 


nt 


‘2 


Time inlorvat, IAgI lo SRG liansler pube t 


0 


ne 


‘3 


Pulso Jiiialion. lAG high 


350 


ne 


«4 


Pulse diiialion. SRG liansloi pulee Ngh 


350 


ne 


'5 


Tmie Hilorval, lAGl lo SRG liansloi (xibo I 


350 


ne 


'6 


Tune inlorval, SRG liaiistor fMitse 1 to SRQ clock pulso T 


70 


ne 


Capacitive load 


OUT pin 


12 




0|ioi.itiMg Iruu air lentparaloie. T/y 


-to 45 


•c 



t TIh) olgiitfak convonltoii, in wlittlt ll>o loabi pobilivo (most noy/itivo) vithiu b niininitKtt. b nso<l in Him Jala tlraul foi clock voH^iye 

kjvols 

^ A'lfiiblmiml IS 101(0110(1 lof Of)lniiul poiloMtiaiicu 



, Texas 
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electrical characteristics over recommended operating range of supply voltage, = - 10X to 40X 



PARAMETER 


MtN TYP< MAX 


UNIT 


Dynamic range (see Mole 2) 


Anliblooming disabled (see Note 3) 


60 


cJB 


AnliWooining enabled 


57 


Ctierge conversbn fed or 


40 


nV/e 


Charge translor effidoncy (eee Mole 4) 


0 00090 0 00096 




Signal response dolsy time, r (see Note 5 end Figuie 5) 


25 


ns 


Gamma (see Note 6) 


0 07 0 06 0 00 




Oirtput teaielarKa 


700 600 


II 


TToise voHage 


\l\ noise (5 kH 2 ) 


370 


tiV//M? 


RerYdom rroise, 1 - tOO Wtf 


70 


Noise equivalent signal 


150 


electrons 


Rejection ratio al 7 l6MHj 


From ADB to output (soe Note 7) 


to 


dB 


F rom SRO to output (soe Note 0) 


37 


Si^jply current. Iqq 


s to 


mA 


CnpacHance 


Image area gete 






SorinI rogislor gate 


25 


Anliblooming gete 


760 



t All lypicel vBlues ete a1 T/y - 25X 



NOT eS' 2 Dynnnik tnngo n - 20 linios Min lo^nrilhm ot Mio nionn noise sigtiol dtvklorl tiy llie snluinlron owtfnif signni 

3 For Mils lesi, Ihe enlQ rloomltiQ gnle mii<;l l>n tiinsnd nl Mio Inloniicdinln lovni 

4 Chmge Irnnslet effirlency h one mliHts Hie fJinipe loss pe^ Menslof In ll»o oolpirt teglslei Tim losi It fKjiloiniod hi Ihe dn»k ueing an 
elorliicnl Inpirt signal 

5 Signal rnspoitSB dolny Inne le Bio lirno Imlwoon llm l.iWhig o<lgn ol llm SFIO dork ptilse and llie oiil|ktl signal vnlH slalo 

6 Ganinin (y) Is R«e vnliio ot Hie erpoftoni In Mm e<^inlrnn below tor Iwe | minis on Hie kiieei poiMon of M»e liansloi fiairtion cirrve (Hits 
value lepiesenls points nae> enimntion) 

( Erposiwe { 2 ) V _ / ChHpul signal ( 2 ) 

Ejiposme (I) I \ Oiilpirt signni (ij 

7 AD0 injoclion iplio Is - 20 timos Ibo log'iiilliin ol Mio oc amfiMiHlo nl Ihe otit|Mil rlivHod by Hio nr. mnpliludo al ADB 

0 SnO rejncHon lallo Is -20 tirnos Ihe Ingniilhm ol Ihe nr, nmfJIliido al Iho onl|>iil divided hy Iho ar nnif«lilnde nl SHO 



, TtXAS 
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Optical characteristics, = 25“C (unless otherwise noted) 



PARAMETER 


MIN TYP MA)( 


UNIT 


No in fitlei 


Mudsiriud al Vy (suo Nolo 10) 


260 


mV/lx 


Se..s.K..,M.Noleg) WHhlRHI*. 


33 


SatiJialion signcii (see Note 1 1) 


AnlAjloomlng drsablod 


400 600 


nrV 


Anlibloomtng onof >lod 


3S0 450 


Bloonimg ovedoad ratio (soe Note 1 2) 


Slrtiie 


5 




Shidtoiod kgtrl 


100 


Output signal nonunMormity ft/2 seturation) (see Note 13) 


10% 20% 




Image aroa well capacity 


150x10^ 


okxdrons 


Oeik current |T/^-2I*C 


0 027 


nA/cm^ 


Dark signal (see Nole 14) 


10 15 


nrV 


Dark signal nonuntformtly for eniKe bold (see Nole 1 5) 


4 1$ 


mV 


Mixiulation lianslet kirKbon 


Hortionlal 


50% 




Vorlkal 


70% 



NOTES: 0 SonsHiviry Is measured el sn inlegialton llmo of 16 667 nis end e sonico lunipoiatitio of 26S6 K A CM SOO Moi (« used 

10 V^j is Iho out(Mi1 volljge Uial lopfesonts thu ltnobliol<l of o|totijtioii of onMioutning V() 1/2 s<jlor;tlion si^itel 

1 1 Saluretion Is 6re coodHion In wtiicti lirrihoi ktcieeso In eif^osiiio d<»os ik> 1 toad to fuiHior incrutiso fn oirt|Mi( srgoal 

12 Blomniug ovuitoed rjbo b lliu ruho of t>looming eiposuio to boliir.tlicNt o>poMuo 

13 Oti1|Kit signal nonnintuiniilyb the lalio of Iho n>silniMm|>tiul lo | tlauldiNoiunro In output sigrtul to thon)oanoiil|Ni1 sigiuW for exposme 
B'tf istud to give \f2 ttiu saluiiition output sign-it 

14 Dork signal lovel is inoasure<t fiom ttie dumniy piiuls 

1 5 Dotk signol nonunHormity w IIk) nianmum pivul to pliol diNuroucu in e d.«k coiKkliou 



PARAMETER MEASUREMENT INFORMATION 




Figure 3. Typical Clock Waveform for lAG and ABG 

Slew RATE BETWEEN 10% end 00% - 70 lo l20V/|i8.t, - 150 ns. I| - 00 i.s 



, Texas 
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PARAMETER MEASUREMENT INFORMATION 




Figure 4 . Typical Clock Waveform for SRG 




Figure 5. SRO end CCD Output Waveforms 
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TYPICAL CHARACTERISTICS 



HORIZONTAL MODULATION 
TRANSFER FUNCTION 
(BARS PARAt LEL TO SERIAL REGISTER) 




Noimalltfld Sptll*l Fftquancy 

—I I I I 

13 I3B Ilf 350 

Spciltl Fr«qutncy - Cycl«i/rmn 



VERPCM. MODULATION 

transfer function 

(BARS PERPENDICULAR TO SERIAL REGISTER) 
1 

































X - 40( 

Voo - 
Ta - 2 
1 


3 to 700 
12 V 
6»C 


nm Mont 


>cKfomatl 


k: tiehl _ 



NormaHtad Spallal Fraquancy 
J I I l_ 



7 3 Ml 311 39 1 

Sptllal Frtquaney - Cyelaa/mm 



FIguro 6 

NOISE SPECTRUM OF OUTPUT AMPl IFIER 



9 

I 




FIgurt 7 

CCD SPECTRAL RESPONSIVITY 




Figure fl 



Incident Wevelengih - nm 
Figure 9 
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Vss — < 

V|A(lnl> 



VOB 

I _ 



P-Oilvw 

TMS3473 






10 

20 



9 Driver 
9N2B840 



TYPICAL APPLICATION DATA 

VaBI 
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• Wv 



YAB(lnl) 



Vab- 



vcc 




User Defined 
Clock OeneiBlor 



Oscltlslor 
J 14 J 



TTTTT 



t)) 



CBIK 

C9YNC 

ClAMP 



SUPPOR! CIRCUITS 


DEVICE 


PACKAGE 


APPLICAIION 


FUNCTION 


SN2B048 


20 pif» NF 


Sottnl d*K/« 


Dirvof lor SRG 


TMS3473 


20 r>*n HF 


Pnrnllel diivm 


Drivw bf 1AO. ABO 


TU503 


16 pin D 


Snniple nod lK>ld 


Thfoo-dtannel sample and hoW 1C 



Figure 10. Typical Application Circuit Diagram 
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MECHANICAL DATA 



r 

^ lie package (or the TC2 1 1 consists ol a ceramic base, glass wirxiow, arxi a 6 lead Irame. The glass window 
is sealed to the package by an epoxy adhesive The package leads are contigiired in a dual In line organization 
and (it into moiinling holes with 2,54 mm (0 1 1nch) center to center spacings 




7,54 (0 297) 
7,14 (6 36 ij 

I® 

I® 

I® 




f'K)1ES A DinitHtsioKS oie IM inHIitnulots 8ikJ |><iiuiiIIioIk ully m kk hus Singlo dtinetisioiin iiio iMintiiitil 
B Tlio coitlei ol Iiio |idrk j<ju mvl Ihu cunlui of Uio imoQO moa h>u moI coiik 

C ThocAst.iof otuMM llio top olilioyl.ibb lollio im.iyo bonsoi smi.tco t$ lypn nlly 1 imii (0 040 iiicli) TImj ykist is lyiitcuily 0 030 Nidi lltkk 
ami hns on jiidoi ol lufiodion ol I 52 
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APPENDIX B MTF PROGRAM 



/*MTF.c find ro and f from the Isf using FFT to the mTF*******/ 

#define TITLE "Measure ro, plot the MTF" 

#define AUTHOR "By W. J. Rail" 

^***********Qyj0^ 0**********************************^ 

/*#include <alloc.h> */ 

#include <lynxx.h> 

#define True 1 
#define False 0 
#define false 0 

#include <stdio.h> 

#include <stdlib.h> 

#include <graph.h> 

#include <math.h> 



#include <sys\types.h> 

,/*int frame_ptr[FRAME_BYTES/2];*/ 
int size=60, xstart=0, ystart=0; 
float lsfx[1 65],lsfmax=0.0,lsfy[1 92]; 

int sign=1, m=8, mm=2, mm2;/*mm is # of points in FFT, mm2=mm/2***/ 
float mtfx[1024],mtfmax=0.0,mtfy[1024], re[1024],im[1024],px[512],py[512j; 
float rocal = 1.0,rox,roy; 

void init(): 

void subexpose(long tm,int shutter); 
void expose(long tm); 
void long_delay(long tm) ; 
void calcjsf (unsigned int far *sub_frame); 
void backgroundQ; 

void fft(int m, int sign, float re[], float im[j); 
void calc_ro(); 
void init_mtf(); 
void plot_mtf(); 

/*to see Isf’s, use these two functions instead of mtf ones*/ 
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void init_graph(); 
void plot_graph(); 

j* main routine */ 

void main() 

{ 

int i; 

for (i=0; i<m-1; i++) mm = 2*mm; 
initJynxxQ; 
initO: 

cooler_off(); 

_setvideomode( _TEXTC80); 

} 

/* Initialize input info */ 

void init() 

{ 

int cool = 0,shutter=0; 
long tm; 

printf("Expose time in ms.:"); 
scant ("%ld",&tm): 
printf("Cooler 0=OFF 1=ON;"); 
scanf("%d",&cool): 
if (cool) 

cooler_on(): 

printf("Use shutter ? (1=Yes/0=No)"); 
scanf("%d", &shutter); 
printf ("Subframe size? "); 
scant ("%d", &size); 

/* printf("\nHit any key to quit:");*/ 

subexpose(tm, shutter) ; 

} 

/* subExpose CCD */ 

void subexpose(long tmjnt shutter) 

{ 

unsigned int *sub_frame; 
int sub_frame_bytes; 

float pxO.pyO; 
int i,j,x,y; 

for (i=0; i<m-1; i++) mm = 2*mm; 
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mm2 = (int) ((float) (mm)/2.0) ; 



sub_frame_bytes = size*size * sizeof(unsigned int); 
sub_frame = (unsigned int*)malloc(sub_frame_bytes); 



if ('shutter) open_shutter(); 



while ( IkbhitO ) 

{ 

if (shutter) 

{ 

expose(tm): 

} 

else 

{ 

cl rood 0: /* if not using shutter 
long_delay(tm); 

} 

x_param = (unsigned char) xstart; 
y_param = (unsigned char) ystart; 
size_param = (unsigned char) size; 
ptr_param = sub_frame; 
digitize_sub_frame(); 
calc_lsf(sub_frame); 
background 0; 






/ 



for (x=xstart+ size-1 ; x>=xstart; x--) re[x] = lsfx[x]; 
fft(m, sign, re, im); 

for (x=0; x<mm; X-1--1-) mtfx[x] = re[x]; 

/*find the x power spectrum***************/ 
pxO = re[0]*re[0] - 1 - im[0]*im[0]; 
for (j=0; j<mm2; j-t--i-) 

{ 

if(mtfx[j] != 0.0 ) if(pxO != 0.0) 

px[j]=sqrt((re[j]*reO]-i-imO]*im[j])/pxO)/mtfx[j]; 

} 

/*reinitialize the FFT arguements before finding the Y FFT*/ 
for (i=0; i<1024; i-i--i-) 

{ 

re[i] = 0.0; 
im[i] = 0.0; 

} 
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for (y=ystart+ size-1; y>=ystart; y-) re[y] = lsfy[y]; 
sign, re, im); 

for (y=0; y<mm; y-f-i-) mtfy[y] = re[y]; 

/*find the y power spectrum***************/ 
pyO = re[0]*re[0] + im[0]*im[0]; 
for (j=0; j<mm2; 

{ 

if(mtfyO] != 0.0) if(py0 != 0.0) 

py[j]=sqrt((reG]*re[j]-l-im[J]*imO])/pyO)/mtfyO]; 

} 

calc_ro(); 

init_mtf(); 

plot_mtf(): 

/*init_graph(): 

plot_graph(); ********* use these to see Isf’s */ 

} 

free(sub_frame); 
cl ose_sh utter 0; 

} 

/* Expose CCD */ 

void expose (long tm) 

{ 

cIrccdO; 

open_shutter(): 

long_delay(tm); 

close_shutterO; 

} 

/* long delay */ 

void long_delay(long tm) 

{ 

delay_param = 1000; 
while (tm > 1000) 

{ 

delay_ms(); 
tm = tm -1000; 

} 

delay_param = (unsigned short)tm; 
delay_ms(); 

} 

/* calculate the line spread functions */ 

void calc_lsf(unsigned int far *sub_frame) 

{ 
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unsigned int pixel; 
int x,y; 

float lsfxmax=0.0,lsfymax=0.0: 
for (x=0; x<165; x++) lsfx[x] = 0.0; 
for (y=0; y<165; y++) lsfy[y]=0.0; 
for (x=size-1+xstart; x>=xstart; x--) 

{ 

for (y = size-1 +ystart; y>=ystart; y--) 

{ 

pixel=*(sub_frame-i- -I-); 
lsfx[x] = lsfx[x] -I- (float)pixel; 
lsfy[y] = lsfy[y] + (float)pixel; 
if(lsfx[x] > Isfxmax) lsfxmax=lsfx[x]; 
if(lsfy[y] > Isfymax) lsfymax=lsfy[y]; 

} 

} 

lsfmax= Isfxmax; 

if(lsfxmax < Isfymax) Isfmax = Isfymax; 

} 

/* Subtract off the background */ 

void backgroundQ 

{ 

int x,y,edge; 

float sumx=0.0,sumy=0.0,ave_x1 ,ave_x2,ave_y1 ,ave_y2,aveback; 
edge=(int)(size/10.0); /*this uses 20% to find background*/ 

/* — Find the average background */ 

for (x=xstart; x<xstart-i-edge; X- 1 -- 1 -) sumx=sumx-i-lsfx[x]; 
ave_x1 =sumx/edge; sumx=0.0; 

for (x=xstart-i- size-edge-1 ; x<xstart-i-size; X- 1 -- 1 -) sumx=sumx-i-lsfx[x]; 
ave_x2=sumx/edge; 

for (y=ystart; y<ystart-i-edge; y-i--i-) sumy=sumy-i-lsfy[y]; 
ave_y1 =sumy/edge; sumy=0.0; 

for (y=ystart-i-size-edge-1 ; y<ystart-i-size; y-i--i-) sumy=sumy-i-lsfy[y]; 
ave_y2=sumy/edge; 

aveback =((ave_x1 -i-ave_x2)-i-(ave_y1 -i-ave_y2))/4.0; 

/*— -subtract off the background */ 

Isfmax = Isfmax - aveback; 

for (x=xstart-i-size-1 ; x>=xstart; x-) 

lsfx[x] = Isfx[x]-((x-xstart)*(ave_x2-ave_x1 )/(size-edge) - 1 - ave_x1 ) ; 
for (y = y start -I- size- 1 ; y>=ystart; y--) 

Isfy [y] = Isfy [y] -((y-ystart) * (ave_y2-ave_y 1 )/(size-edge) - 1 - ave_y 1 ) ; 

/* printf(”\nThe average background is %An",aveback);*/ 
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} 



void fft(int m, int sign, float re[], float im[]) 

{ 

int n, j, j1, ndiv2, n1, n2, k, i, ip, npts; 
int le, leO, lei , I; 

float pts, t, ure, uim, wre, wim, tre, tim; 
double ang, pi; 

pi = 4,0*atan(1.0); 

n = (int)(pow(2.0, (double)m) + 1.0e-10); 

j = 1; 
j1 = 0; 
n1 = n - 1 ; 
n2= n - 2; 
ndiv2 = n/(int)2; 



for(i=0; i<=n2; i++) 

{ 

if(i < j1) 

{ 

t = reyi]; /* you might do better here with pointer operations 

re[j1] = re[i]; 
re[i] = t; 
t = im[j1]: 
imyi] = im[i]: 
im[i] = t; 

} 

k = ndiv2; 
while(k < j) 

{ 

j -= k; 
k /= 2; 

} 

j + = k; 



le — 1 ; 

for(l=1; l<=m; I++) 

{ 
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ure = (float) 1; 
uim = (float)O; 
ang = pi/(double)le; 
wre = (float)cos(ang); 
wim = (float)sin(ang); 
Ie0=le; 

Ie1=le-1; 
le +=le; 



forO=0; j<=le1; j + +) 

{ 

for(i=j; i<=n1; i+=le) 

{ 

ip = i + leO; 

tre = re[ip]*ure - im[ip]*uim; 
tim = re[ip]*uim + im[ip]*ure; 
re[ip] = re[i] - tre; 
im[ip] = im[i] - tim; 
re[i] += tre; 
im[i] += tim; 

} 

t = ure*wre - uim*wim; 
uim = ure*wim + uim*wre; 
ure = t; 



} 

} 

if(sign<0) 

{ 

pts = 1/(float)n; 

for(i=0; i<n; i++) 

{ 

re[i] *= pts; 
im[i] *= pts; 

} 

} 



} 



void calc_ro() 

{ 

int j; 

float z,zold,pj, slope; 



if (px[1] > 1.0) if (px[1] <= 0.0) printf(W DATA POINT BAD"); 
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else 

{ 

z = log (-log(px[1])): 
for (j=2; j<mm2; j ++) 

{ 

pj = px[j]: 

zold = z; 
if (pj<1.0) 

{ 

z = log(-log(pj)): 

slope = (z-zold)/log((float)(j)/(float)0-1)): 
if (z > 0.0) goto ro; 

} 

} 

} 

ro: rox = rocal*exp(-zold/slope + log(j-1)); 

if (py[1] > 1.0) if (py[1] <= 0.0) printf("\nY DATA POINT BAD") 
else 
{ 

z = log (-log(py[1])): 
for (j=2; j<mm2; j++) 

{ 

pj = py[j]; 

zold = z; 
if (pj<1.0) 

{ 

z = log(-log(pj)); 

slope = (z-zold)/log((float)(j)/(float)(j-1)); 
if (z > 0.0) goto ro2; 

} 

} 

} 

ro2; roy = rocal*exp(-zold/slope + log(j-1)); 

} 

/* initialize for plotting-MTF’s */ 

void init_mtf() 

{ 

int TRUE=1; 
int x,y: 

for (x=0; x<mm; x++) 
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if(mtfx[x] > mtfmax) mtfmax=mtfx[x]; 
for (y=0; y<mm; y++) 

if(mtfy[y] > mtfmax) mtfmax=mtfy[y]; 

setvideomode( _VRES16C0L0R); 
setviewport(0, 39, 639,439); 

_setwindow(TRUE,0,0,mm+mm, mtfmax); 
setcolor(15); 
moveto_w(0.,0.); 



/* plot modulation transfer functions */ 

void plot_mtf() 

{ 

int x,y; 

for (x=mm2; x<mm; x+ + ) 

{ 

_lineto_w( (doubi e) (x-mm2) , (doubi e) mtfx[x]) ; 

} 

for (x=0; x<mm2; x++) 

{ 

_lineto_w((double)(x+mm2),(double)mtfx[x]); 

} 

for (y=mm2; y<mm; y++) 

{ 

_lineto_w((double)(mm2+ y),(double)mtfy[y]); 

} 

for (y=0; y<mm2; y++) 

{ 

Jineto_w((double)(mm+ mm2 + y),(double)mtfy[y]); 

} 

_settextposition(1 ,1 ); 
printf(”%s\n %s",TITLE,_DATE_); 

_settextposition (4,1); 
printfC'rox = %f',rox); 

_settextposition(4,38); 
printf("roy = %f',roy); 

_settextposition(40,30); 

_outtext(”Hit any key to exit:"); 
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/* getchO; there is also a getch in subexpose, only use one*/ 

} 



/* initialize for plotting-Isf’s */ 

void init_graph() 

{ 

int TRUE=1; 

_setvideomode( _VRES16COLOR); 

_setviev\/port(0,39,639,439); 

_setwindow(TRU E, 0, 0, size + size, Isfmax) ; 

_setcolor(15); 

_moveto_w(0.,0.); 

} 

/* plot line spread functions */ 

void plot_graph() 

{ 

int x,y; 

for (x=xstart; x<xstart+size; x++) 

{ 

_lineto_w((double)(x-xstart),(double)lsfx[x]): 

} 

for (y=ystart; y<ystart+size; y++) 

{ 

_lineto_\A/((double)(size+ y-ystart),(double)lsfy[y]); 

} 

_settextposition(1 ,1 ); 
printf ("%s\n %s\n\n",TITLE,_DATE_) ; 

_settextposition (40, 30) ; 

_outtext("Hit any key to exit:"); 

/* getchO: there is also a getch in subexpose, only use one*/ 

} 
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APPENDIX C ARTIFICIAL STAR PROGRAM 



/*Uses an artificial star to test the MTF.c program*has Isf plot also**/ 
#define TITLE "Plot an artificial star MTF" 



#define AUTHOR 

/********Quick C 



"By W. J. Rail" 

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 



#include <stdio.h> 
#include <stdlib.h> 
#include <graph.h> 
#include <math.h> 



#include <sys\types.h> 



int xsize=165, ysize=192, xstart=0, ystart=0; 
float lsfx[1 65],lsfmax=0.0,lsfy[1 92]; 

int sign=1, m=8, mm=2, mm2; 

float mtfx[1 024],mtfmax=0.0,mtfy[1 024], re[1 024],im[1 024],px[51 2],py[51 2]; 
float rocal = 1.0,rox,roy; 



void make_star(); 
void background 0; 

void fft(int m, int sign, float re[], float im[]); 
void calc_ro(); 
void init_mtf(); 
void plot_mtf(); 

/*to see Isf’s use these two functions instead of mtf ones*/ 
void init_graph(); 
void plot_graph(); 



/* rnain routine */ 

void main() 

{ 

float pxO,pyO; 
int i,j,x,y; 

for (i=0; i<m-1; i++) mm = 2*mm; 
mm2 = (int) ((float) (mm)/2.0); 
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make_star(); 
background 0; 

for (x=xstart+xsize-1; x>=xstart; x--) re[x] = lsfx[x]; 
fft(m, sign, re, im); 

for (x=0; x<mm; x++) mtfx[x] = re[x]; 

/*find the x power spectrum***************/ 
pxO = re[0]*re[0] + im[0]*im[0]; 
for (j=0; j<mm2; j ++) 

{ 

if(mtfx[j] != 0.0 ) if(px0 != 0.0) 

pxO ] = sqrt((re [j] *re [j] + imO] *im[j])/px0)/mtfx[j] : 

} 

/*reinitialize the FFT arguements before finding the Y FFT*/ 
for (i=0; i<1024; i + +) 

{ 

re[i] = 0.0; 
im[i] = 0.0; 

} 

for (y=ystart+ysize-1 ; y>=ystart; y--) re[y] = lsfy[y]; 
fft(m, sign, re, im); 

for (y=0; y<mm; y++) mtfy[y] = re[y]; 

/*find the y power spectrum***************/ 
pyO = re[0]*re[0] + im[0]*im[0]; 
for (j=0; j<mm2; j ++) 

{ 

if(mtfyO] != 0.0) if(py0 != 0.0) 
pyO]=sqrt((re[j]*reO]+im[j]*im|j])/pyO)/mtfy[j]; 

} 

calc_ro(); 

init_mtf(); 

plot_mtf(); 

/* init_graph(); 
plot_graph(); */ 

_setvideomode( _TEXTC80); 

} 

/* make star */ 

void make_star() 



64 



int x,y: 

float pixel, lsfxmax=0.0,lsfymax=0.0; 
for (x=xstart+xsize-1 : x>=xstart; x--) 

{ 

for (y=ystart+ysize-1 ; y>=ystart; y--) 

{ 

pixel = 4095*exp (-0. 1 * (float) (((x-82) * (x-82) + (y-96) * (y-96)) ) ) ; 

lsfx[x] = lsfx[x] + pixel; 

lsfy[y] = lsfy[y] + pixel; 

if(lsfx[x] > Isfxmax) lsfxmax=lsfx[x]; 

if(lsfy[y] > Isfymax) lsfymax=lsfy[y]; 

} 

} 

I sfm ax = Isfxmax; 

if(lsfxmax < Isfymax) Isfmax = Isfymax; 

} 

/* Subtract off the background */ 

void background 0 

{ 

int x,y,edge; 

float sumx=0.0,sumy=0.0,ave_x1,ave_x2,ave_y1,ave_y2,aveback; 
edge=(int)(ysize/10.0); /*this uses 20% to find background*/ 

/* — Find the average background */ 

for (x=xstart; x<xstart+edge; x++) sumx=sumx+lsfx[x]; 
ave_x1 =sumx/edge; sumx=0.0; 

for(x=xstart+xsize-edge-1; x<xstart+xsize;x++)sumx=sumx+lsfx[x]; 
ave_x2=sumx/edge; 

for (y=ystart; y<ystart+edge; y++) sumy=sumy+lsfy[y]; 
ave_y1 =sumy/edge; sumy=0.0; 

for (y = ystart + ysize-edge-1 ; y < ystart + ysize ; y+ + ) 
sumy=sumy+lsfy[y]; 
ave_y2=sumy/edge; 

aveback = ((ave_x1 +ave_x2) +(ave_y1 +ave_y2))/4.0; 

/*— -subtract off the background */ 

Isfmax = Isfmax - aveback; 

for (x=xstart-i-xsize-1 ; x>=xstart; x--) 

lsfx[x] = Isfx[x]-((x-xstart)*(ave_x2-ave_x1 )/(xsize-edge) -f ave_x1 ) ; 
for (y=ystart+ysize-1 ; y>=ystart; y-) 

Isfy [y] = Isfy [y]-((y-ystart)*(ave_y2-ave_y 1 )/(ysize-edge) -i-ave_y 1 ); 

/* printf("\nThe average background is %An", aveback);*/ 

} 
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th© FFX *******^********************** 
void fft(int m, int sign, float re[], float im[]) 

{ 

int n, j, j1, ndiv2, n1, n2, k, i, ip, npts; 
int le, leO, lei, I; 

float pts, t, ure, uim, wre, wim, tre, tim; 
double ang, pi; 

pi = 4.0*atan(1.0); 

n = (int)(pow(2.0, (double)m)+1 .Oe-10); 

j = 1: 
j1 = 0; 
n1 = n - 1 ; 
n2= n - 2; 
ndiv2 = n/(int)2; 



for(i=0; i<=n2; i++) 

{ 

if(i < j1) 

{ 

t = reyi]; /* you might do better here with pointer operations 

re[j1] = re[i]; 
re[i] = t; 
t = im[j1]: 
imOl] = im[i]; 
im[i] = t; 

} 

k = ndiv2; 
while(k < j) 

{ 

j -= k; 
k/= 2; 

} 

j + = k; 



le = 1; 

for(l=1; l<=m; I + +) 

{ 

ure = (float) 1; 
uim = (float)O; 
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ang = pi/(double)le; 
wre = (float)cos(ang); 
wim = (float)sin(ang): 
Ie0=le; 

Ie1 =le-1 : 
le +=le; 



for(j=0; j<=le1; j++) 

{ 

for(i=j; i<=n1; i+=le) 

{ 

ip = i + leO; 

tre = re[ip]*ure - im[ip]*uim; 
tim = re[ip]*uim + im[ip]*ure; 
re [ip] = re[i] - tre; 
im[ip] = im[i] - tim; 
re[i] += tre; 
im[i] += tim; 

} 

t = ure*wre - uim*wim; 
uim = ure*wim + uim*wre; 
ure = t; 



} 

} 

if(sign<0) 

{ 

pts = 1/(float)n; 

for(i=0; i<n; i + +) 

{ 

re[i] *= pts; 
im[i] *= pts; 

} 

} 

} 

i|ofp rp**** ****★**★★****★**★★ *★★*★★**** 



void calc_ro() 

{ 

int j; 

float z,zold,pj, slope; 



if (px[1j > 1.0) if (px[1j <= 0.0) printf("\nX DATA POINT BAD"); 
else 
{ 



67 



z = log {-log(px[1])); 
for (j=2; j<mm2; j ++) 

{ 

pj = pxD]: 
zold = z; 
if (pj<1.0) 

{ 

z = log(-log(pj)); 

slope = (z-zold)/log((float)(i)/(float)(j-1)): 
if (z > 0.0) goto ro; 

} 

} 

} 

ro: rox = rocal*exp(-zold/slope + log(j-1)); 

if (py[1] > 1.0) if (py[1] <= 0.0) printf("\nY DATA POINT BAD") 
else 
{ 

z = log (-log(py[1])); 
for (j=2; j<mm2; j + +) 

{ 

PJ = pyQ]: 
zold = z; 
if (pj<1.0) 

{ 

z = log(-log(pj)); 

slope = (z-zold)/log((float)(j)/(float)(j-1)); 
if (z > 0.0) goto ro2; 

} 

} 

} 

ro2; roy = rocal*exp(-zold/slope + log(j-1)); 

} 

/* initialize for plotting--MTF’s *! 

void init_mtf() 

{ 

int TRUE=1; 
int x.y: 

for (x=0; x<mm; x++) 

if(mtfx[x] > mtfmax) mtfmax=mtfx[x]; 
for (y=0; y<mm; y++) 
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if(mtfy[y] > mtfmax) mtfmax=mtfy[y]; 

setvideomode{ _VRES16C0L0R); 
setviewport(0, 39, 639, 439) ; 

_setwindow(TRUE,0,0,mm+mm, mtfmax); 
setcolor(15): 
moveto_w(0.,0.): 



/* plot modulation transfer functions */ 

void plot_mtf() 

{ 

int x,y; 

for (x=mm2; x<mm; x++) 

{ 

_lineto_w((double)(x-mm2),{double)mtfx[x]); 

} 

for (x=0; x<mm2; x++) 

{ 

_lineto_w({double)(x+mm2),(double)mtfx[x]): 

} 

for (y=mm2; y<mm; y++) 

{ 

_lineto_w{{double)(mm2+ y),(double)mtfy[y]); 

} 

for (y=0; y<mm2; y++) 

{ 

_lineto_w{(double)(mm+ mm2 + y),{double)mtfy[y]); 

} 



_settextposition(1 ,1): 
printf{”%s\n %s",TITLE,_DATE_); 

_settextposition (4,1); 
printfC'rox = %f,rox); 
_settextposition{4,38); 
printfC'roy = %f,roy); 

settextposition (40, 30) ; 
outtext("Hit any key to exit;"); 
getchO; 
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} 



/* initialize for plotting — LSF’s — */ 

void init_graph() 

{ 

int TRUE=1; 

_setvideomode( _VRES16COLOR); 
_setviewport(0,39,639,439); 

_setwindow(TRUE,0,0,xsize+ysize,lsfmax): 

_setcolor(15): 

_moveto_w(0.,0.): 

} 

/* plot line spread functions */ 

void plot_graph() 

{ 

int x,y; 

for (x=xstart; x<xstart+xsize; x++) 

{ 

Jineto_w((double)(x-xstart),(double)lsfx[x]): 

} 

for (y=ystart; y<ystart+ysize; y++) 

{ 

_lineto_v\/((double) (xsize + y-ystart), (double)lsfy [y]) ; 

} 

_settextposition(1 ,1): 
printf("%s\n %s\n\n", TITLE, _DATE_); 

_settextpositi on (40, 30) ; 

_outtext("Hit any key to exit:"); 
getchO: 

} 
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APPENDIX D JITTER PROGRAM 



/*jitter.c**********finds the cetroid Jitter, ro and f, plots isf******/ 
#define TITLE "Jitter calculations for ro and f, plots Isf’s" 



#define AUTHOR "By W. J. Rail" 



/*#include <alloc.h> */ 
#include <lynxx.h> 



#define True 1 
#define False 0 
#define false 0 



#include <stdio.h> 
#include <stdlib.h> 
#include <graph.h> 
#include <math.h> 
#include <time.h> 



#include <sys\types.h> 

/*int frame_ptr[FRAME_BYTES/2];*/ 
int size=160, xstart=0, ystart=0; 
float lsfx[165],lsfmax=0.0,lsfy[192]; 

float xcsum =0.0, xc2sum=0.0, ycsum=0.0, yc2sum=0.0; /*for jitter calcs*/ 

int n=1; /*number of centroids between each jitter calc*/ 

float dia=1.0, focal = 1000.0, k=12.57E+6;/*wave number, 500nm*/ 

float xjit,yjit, rocal = 1.0, rox,roy,fx,fy; 

void init(); 

void subexpose(long tm,int shutter); 
void expose (long tm); 
void long_delay(long tm) ; 
void calcjsf (unsigned int far *sub_frame); 
void backgroundO; 
void centroidO: 
void jitO: 
void jit_rof(); 
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void init_graph(); 
void plot_graph{); 

— main routine— 



/* 

void main() 

{ 

initJynxxO: 
initO; 

cooler_off(); 

_setvideomode( _TEXTC80); 



*/ 



/* Initialize program info */ 

void initQ 

{ 

int cool = 0,shutter=0; 
long tm; 

printf ("Expose time in ms.:"); 
scant ("%ld",&tm): 

printf("Cooler 0=OFF 1=ON:"); 
scanf("%d",&cool): 
if (cool) cooler_on(); 

printf("Use shutter ? (1 =Yes/0=No)"); 
scanf("%d", &shutter); 

printf("Enter size of subframe in pixels:"); 
scanf("%d",&size); 

printf("\nHit any key to quit:"); 



subexpose(tm, shutter); 

} 

/* subExpose CCD */ 

void subexpose(long tmjnt shutter) 

{ 

unsigned int *sub_frame; 
int sub_frame_bytes, i; 
double deltime,time1 ,time2; 
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sub_frame_bytes = size*size * sizeof(unsigned int); 
sub_frame = (unsigned int*)malloc(sub_frame_bytes); 

if (Ishutter) open_shutter(); 

while ( IkbhitO ) 

{ 

for(i=0; i<n; i++) 

{ 

if (shutter) expose (tm); 
else 
{ 

clrccdQ; /*if not using shutter*/ 
long_delay(tm): 

} 

x_param = (unsigned char) xstart; 
y_param = (unsigned char) ystart; 
size_param = (unsigned char) size; 
ptr_param = subjrame; 
digitize_sub_frame(): 

calc_lsf (subjrame); 

backgroundQ; 

centroidO; 

} 

jit(): 

jit_rof(); 

init_graph(); 

plot_graph(); 

} 

free(subjrame); 

close_shutterO; 

} 

/* Expose CCD */ 

void expose(long tm) 

{ 

clrccdQ; 
open_shutter(); 
long_delay(tm); 
cl ose_sh utter 0; 
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/*-- long delay */ 

void long_delay(long tm) 

{ 

delay_param = 1000; 
while (tm > 1000) 

{ 

delay_ms(); 
tm = tm -1000; 

} 

delay_param = (unsigned short)tm; 
delay_ms(); 



/* calculate the line spread functions */ 

void calc_lsf(unsigned int far *sub_frame) 

{ 

unsigned int pixel; 
int x,y; 

float lsfxmax=0.0,lsfymax=0.0; 
for (x=0; x<165; x+-i-) lsfx[x]=0.0; 
for (y=0; y<165; y++) lsfy[y]=0.0; 
for (x=size-1 +xstart; x>=xstart; x-) 

{ 

for (y = size-1 +ystart; y>=ystart; y-) 

{ 

pixel=*(sub_frame++); 
lsfx[x] = lsfx[x] + (float) pixel; 
lsfy[y] = lsfy[y] + (float)pixel; 
if(lsfx[x] > Isfxmax) lsfxmax=lsfx[x]; 
if(lsfy[y] > Isfymax) lsfymax=lsfy[y]; 

} 

} 

lsfmax= Isfxmax; 

if(lsfxmax < Isfymax) Isfmax = Isfymax; 

} 

/* Subtract off the background */ 

void backgroundQ 

{ 

int x,y,edge; 

float sumx=0.0,sumy=0.0,ave_x1 ,ave_x2,ave_y1 ,ave_y2,aveback; 
edge= (int) (size/1 0.0); /*this uses 20% to find background*/ 

/* — pind the average background */ 
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for (x=xstart; x<xstart+edge; x++) sumx=sumx+lsfx[x]; 
ave_x1 =sumx/edge; sumx=0.0; 

for (x=xstart4-size-edge*1 : x<xstart+size; x++) sumx=sumx+lsfx[x]; 
ave_x2=sumx/edge; 

for (y=ystart; y<ystart+edge; y++) sumy=sumy+lsfy[y]; 
ave_y1 =sumy/edge; sumy=0.0; 

for (y=ystart+size-edge-1 ; y<ystart+size; y++) sumy=sumy+lsfy[y]; 
ave_y2=sumy/edge; 

aveback = ((ave_x1 + ave_x2) + (ave_y 1 + ave_y2))/4.0; 

/* — subtract off the background */ 

Isfmax = Isfmax - aveback; 

for (x=xstart+ size-1 ; x>=xstart; x--) 

lsfx[x] = Isfx[x]-((x-xstart)*(ave_x2-ave_x1 )/(size-edge) + ave_x1 ) ; 
for (y=ystart+ size-1; y>=ystart; y--) 

Isfy[y]=lsfy[y]-((y-ystart)*(ave_y2-ave_y1 )/(size-edge)+ave_y1 ); 
r printf^nThe average background is %f\n",aveback);*/ 

} 

/* Find the centroids */ 

void centroidQ 

{ 

int x,y; 

float xcent,ycent,sumlx=0.0,sumly=0.0,sumxlx=0.0,sumyly=0.0; 

for (x=xstart-f size-1 ; x>=xstart; x-) 

{ 

for (y =ystart 4- size-1 ; y>=ystart; y-) 

{ 

sumxlx = sumxlx-f- lsfx[x]*(float)(x); 
sumix = sumix + lsfx[x]; 
sumyly = sumyly+ lsfy[y]*(float)(y); 
sumly = sumly -f- lsfy[y]; 

} 

> 

xcent=sumxlx/sumlx; ycent=sumyly/sumly; 

/*printf("\n Centroid is xcent=%f \n ycent=%f",xcent,ycent);*/ 

/* sum centroids for jitter calcs****************/ 
xcsum = xcsum + xcent; xc2sum = xc2sum -f- xcent*xcent; 
ycsum = ycsum -f- ycent; yc2sum = yc2sum 4 - ycent*ycent; 



void jitO 



jitter 



************************** 



/ 
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{ 

xjit = sqrt((xc2sum-xcsum)/(float)(n)); 
yjit = sqrt((yc2sum-ycsum)/(float)(n)); 
xcsum=0.0; xc2sum=0.0; ycsum=0.0; yc2sum=0.0; 
printf('An Jitter is xjitter=%f yjitter=%f, xjit, yjit); 

} 

/*******************Q 3 lQy|g^g |.Q gpjj f *f|-onn the jitter****************/ 
void jit_rof() 

{ 

float a,b; 

a = xjit * k / focal; 

b = 1.418 * po\A^(dia, 0.333) * pow(a,2); 
rox = rocal / pow(b,0.6); 

a = yjit * k / focal; 

b = 1.418 * pow(dia, 0.333) * pow(a,2); 
roy = rocal / pow(b,0.6); 

fx = xjit; 
fy = yjit*roy; 

} 

/* initialize for plotting */ 

void init_graph() 

{ 

int TRUE=1; 

_setvideomode( _VRES16C0L0R); 

_setviewport(0,39,639,439); 

_setwindo\A^(TRUE,0,0,size+size,lsfmax); 

_setcolor(15); 

_moveto_w(0.,0.); 

} 

/* plot line spread functions */ 

void plot_graph() 

{ 

int x,y; 

for (x=xstart; x<xstart+size; x++) 

{ 

_lineto_w((double)(x-xstart),(double)lsfx[x]); 
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} 

for (y=ystart; y<ystart+size; y++) 

{ 

Jineto_w((double)(size+ y-ystart),(double)lsfy[y]) 

} 

_settextposition(1 ,1); 
printf ("%s\n %s\n\n",TITLE._DATE_) ; 

_settextposition (4,1); 
printfC'rox = %g\nfx = %g",rox,fx); 

_settextposition (4,38) ; 
printf("roy = %g",roy); 

_settextposition (5,38) ; 
printf ("fy = %g",fy); 



_settextposition (40, 30) ; 

_outtext("Hit any key to exit:"); 

/* getchO: subexpose also has one, only use one*/ 

} 
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APPENDIX E FRAME RATE MEASUREMENT CODE 



/*time.c****times digitize, calcisf and background, cetroid, jitter*/ 
#define TITLE “Time calculations" 



#define AUTHOR 

/**Quick C 
#include <time.h> 



"By W. J. Rail" 

******Sample of time code* 






/ 



while ( IkbhitO ) 

{ 

time1 = clockO: /* start clock */ 

for(i=0j i<128j i + +) 

{ 

digitize_sub_frame(): 
calcjsf (sub_f rame) ; 
background 0; 
centroid (); 



} 

jit(i); 



time2=clock(); /* stop clock */ 
deltime = (time2-time1)/CLK_TCK; 

printf("\nTime for 128 frames =%f sec", deltime); 

} 
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